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ABSTRACT 

Time resolved spectra of /3 Cep yield the following average results: Tes =24000 ± 250 K, 
log5 = 3.91 ±0.05 and ^ = 8.1 ±0.9 km s"^ N, O, Ne, Al, Si and S abundances are 
solar while C, Mg and Fe are slightly under-abundant. Pulsational amplitudes of 
ATcff '--^ 700 K and Alogg~0.2 dex are found from H/3. The metal lines give simi- 
lar amplitudes but centred on Tcff '--^ 25000 K. An upper limit of 1.0 km s~^ to the 
variability of the microturbulence is derived from the Si in triplet at 455 nm. The 
radial velocity amplitude derived from the core of H/3 is ~ 15% greater than that from 
the metal lines. 
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1 INTRODUCTION 

The capability of oscillations to probe the internal structure 
of stars is fully testified by solar studies. For this reason, 
pulsating stars, that are spread in the HR diagram, are ex- 
pected to largely improve our knowledge of stellar evolution 
and physical assu mptions in models. 

lOsakil (|l97ll ) has shown that both radial and non-radial 
modes of pulsating stars are observationally distinguishable 
via the analysis of line profile variation (LPV) . From this pi- 
oneering work, there has been an increasing interest in LPV 
modeling driven by the improve ment in resolution and effi- 
ciency of modern spectrographs: iLee fc Said (Il990l') inc luded 
the influence of stellar rotation. Briquet et al.j ( 200 ll ) con- 
sidered photospheric spots, Acrts ct al. (1994) performed a 
mode identification v ia th e moments of the line profile or 
iTelting fc Schriiverd (|l997l ) via the Doppler imaging tech- 
nique. 



iDupret et al.l (|2002l ) noted that a pure kinematic ap- 
proach is not enough for LPV modeling and that a realistic 
representation of the strength and broadening of spectral 
lines were necessary. They considered the atmosphere of a 
pulsating star as a perturbed status of the hydrostatic equi- 
librium and represented the temperature perturbation of at- 
mospheric layers with a series expansion in terms of effective 
temperature and gravity. 

This operative method does not clearly represent the 
physics of /3 Cep stars, pulsation modes do not change the 
energy production and the very small (1%) stellar radius 
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variation IjAerts et al.lll994h cannot modify neither the ef- 
fective temperature nor the surface gravity in an observable 
way. However along a pulsation cycle, the compression and 
expansion of photospheric layers change the dependence of 
temperature and pressure on optical depth and let the effec- 
tive temperature appear variable, with the pulsation period, 
if determined in a limited spectral interval . A phenomenon 
i ndeed observed in some /3 Cephei stars bv lDe Ridder et al.l 
(J2002t ) on the basis of a photometric approach. These au- 
thors derived a Toff variability up to 1000 K with the pul- 
sation period, from Geneva colours. 

With the density and temperature, the variation of 
atomic level populations should also result in apparent sur- 
face gravity variations. In fact, there is an equivalence be- 
tween stellar surface gravity and ele ctron press ure in the 
framework of spectral line formation (|Gravll 19921 ). 

A further contribution to the strength and profile of 
spectral lines comes from the microturbulence (^), that in 
principle could be modified along compression and expan- 
sion phases. 

In this paper we analyse a set of time resolved spectra 
of (3 Cep (= HD 205021 = HR8238), which is generally ac- 
cepted to be the proto-type of a class of hot pulsating stars, 
with the aim of ascertaining the apparent variability of the 
effective temperature, gravity and microturbulence with the 
pulsation period. These variations would strongly constrain 
any LPV method, as the one proposed by Dupret and col- 
laborators. 
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Table 1. For each spectrum we report the Heliocentric Julian 
Date, the radial velocity measured from the Si III triplet at 4550A, 
the Teff and log g determined through the H/3 line. 



HJD 


Vrad 


Toff 


logs 


(2554041.0+) 


(km s-l) 


(K) 




0.2510 


-19.70 ±0.97 


24500 ± 310 


3.95 ±0.04 


0.2587 


-17.86 ±0.85 


24150 ± 240 


3.87 ±0.05 


0.2674 


-16.68 ±0.40 


24300 ± 240 


3.89 ±0.05 


0.2746 


-13.67 ±0.49 


24250 ± 220 


3.90 ±0.05 


0.2819 


-10.01 ±0.35 


24150 ± 220 


3.88 ±0.05 


0.2892 


-6.86 ±0.28 


24050 ± 210 


3.88 ±0.05 


0.2964 


-3.95 ±0.40 


23800 ± 200 


3.81 ±0.05 


0.3040 


0.76 ±0.19 


23700 ± 160 


3.78 ±0.04 


0.3113 


4.89 ±0.30 


23750 ± 170 


3.78 ±0.04 


0.3186 


8.18±0.82 


23700 ± 160 


3.77 ±0.04 


0.3261 


10.63 ±0.55 


23850 ±210 


3.81 ±0.05 


0.3334 


11.97±0.64 


23850 ±210 


3.85 ±0.06 


0.3406 


12.18 ±0.26 


23850 ± 220 


3.84 ±0.06 


0.3515 


10.70 ±0.22 


24150 ±210 


3.88 ±0.05 


0.3588 


8.24 ±0.50 


24400 ± 220 


3.92 ±0.04 


0.3660 


4.95 ±0.28 


24200 ± 210 


3.89 ±0.05 


0.3733 


1.75 ±1.01 


24300 ± 240 


3.90 ±0.05 


0.3806 


-1.00 ±0.99 


24400 ± 260 


3.92 ±0.04 


0.3878 


-5.13 ±1.42 


24250 ± 250 


3.89 ±0.06 


0.3955 


-8.74 ±2.13 


24200 ± 250 


3.86 ±0.06 


0.4028 


-12.14 ±1.04 


24200 ± 250 


3.89 ±0.06 


0.4100 


-15.52±0.71 


24350 ± 250 


3.91 ±0.06 


0.4173 


-18.43 ±0.29 


24350 ± 240 


3.90 ±0.05 


0.4246 


-20.11 ±0.94 


24350 ± 210 


3.92 ±0.04 


0.4318 


-20.13 ±0.88 


24300 ± 210 


3.90 ±0.05 


0.4475 


-19.34 ±1.26 


24050 ± 210 


3.87±0.05 


0.4548 


-15.45 ±2.05 


24050 ± 230 


3.87±0.05 


0.4620 


-14.02 ±1.02 


23850 ±210 


3.83 ±0.05 


0.4693 


-11.86 ±1.57 


23900 ±210 


3.84 ±0.05 


0.4766 


-8.07 ±0.89 


23800 ± 190 


3.81 ±0.05 


0.4838 


-4.53 ±0.92 


23800 ± 190 


3.81 ±0.05 



2 OBSERVATIONS AND DATA REDUCTION 

Spectroscopic observations of 13 Cep have been carried out at 
the 91 cm telescope of the INAF - Osservatorio A strofisico di 
Cata nia at R= 20 000 in the range 4300-6800 A (jCatanzard 
l2008l ). 

Spectra were acquired during the night of November, 
the 1"* 2006, starting at JD = 2454041.2510, for a total of 
31 exposures. Taking into account the pulsational period 
P = 0.19048678 days l|Chapelliej 119851 ) . we set the expo- 
sure time to 10 minutes, « 4% of the pulsational period, 
to avoid phase smearing effects. The stellar spectra, cali- 
brated in wavelength and with the continuum normalised 
to a unity level, were obtained using standard data re- 
duction procedures for sp ectroscopic obse rvations within 
the NOAO/IRAF package (ICatanzardl2008l ). The signal-to- 
noise ratio of our spectra always resulted above 200. 

Instrumental systematic errors on radial velocities have 
been corrected observing stars with constant and well-known 
ra dial velocity t aken from the list of standard stars published 
bv lUdrv et all (|l999l ): HD 186791 and HD 206778. 
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Figure 1. From bottom to top, the radial velocity measured from 
Sim triplet, effective temperature and gravity obtained by fitting 
H/3 line profiles, and nitrogen abundance are plotted as a function 
of time. Vertical dashed lines represent the instants of maximum 
and minimum radius. The over-imposed sinusoidal fit of veloci- 
ties has been cal culated adopting pulsational period derived by 
[ChaEeiiie3 lUii^). 



3 RADIAL VELOCITY VARIABILITY 

Radial velocities (Vrad) of /3 Cep, as it is common for 
this class of pulsating stars, have been measured by 
means of the Sim triplet at 4552, 4567 and 4574 A. These 
are st rong lines almost ins ensitive to temperature varia- 
tions (jDe Ridder et al.ll2002l ) and presenting no significative 
blending and. 

Errors on Vrad have been computed assuming that the 
main source of error is the statistical noise in the observed 
spectrum. Velocities are reported in Tab. [T] together with 
the respective Heliocentric Julian Date and plotted in the 
bottom panel of Fig.[l] The parameters of the over-imposed 
sinusoidal fit have been com puted adopting the pulsational 
period bv lChapellieiJ l|l985r ). The results of fitting procedure 
were: an average value 70 = —4.8 ± 0.1 km s^^ and amplitude 
K =16.3 ±0.2 km s-^ 



4 AVERAGE STELLAR PARAMETERS 

The first step to highlight the variability of effective temper- 
ature (Teff), surface gravity (logg) and microturbulence (^) 
in (3 Cep is the ascertainme nt of their average behaviour. 

iLeone fc Manfrg l|l997l ) have demonstrated that the ef- 
fective temperature and gravity inferred from Balmer lines 
have to be determined simultaneously with abundances. 
These define the metal opacity and with i t the dependence 
of_the temperature on the optical depth (JLeone fc Manfrel 
[l99a). Moreover, the abundances define the Stark broad- 
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ening of Balmer lines through the electron pressure. Simi- 
larly, the effective temperature and gravity of a star have 
to be determined simultaneously with the microturbulence. 
This velocity field defines the metal opacity through the line- 
blocking. This is an effect particularly important for a hot 
star whose flux is peaked in the ultraviolet where the largest 
concentration of a tomic transitions appear. As to /3Cep, 
iMorel et all (| 200d ) determined a rather high (6±3 km s"^) 
microturbulence certainly able to largely desaturate metal 
lines. 

To establish the atmosphere parameters of /3 Cep, 
we implemented the iterative procedure presented in 
iLeone fc Manfrg 11997]) to include the microturbulence. In 
particular we compare the observed and theoretical pro- 
files of H[7 and H/3 lines by minimising the residuals. 
Abundances are derived froma comparison between the ob- 
served average spectrum and SYNTHE spectra. Following 
lErspamer fc NorthI (J2003l ). the observed spectrum has been 
divided to intervals of about 100 A and abundances are 
given by the average values in different intervals. Errors rep- 
resent the standard deviations. The microturbulence veloc- 
ity is determined demanding that the abundances derived 
from 14 N ll and 17 O ll unblended lines are independent of 
equivalent widths. In this case, t he abundances are deri ved 
by means of the WIDTH9 code (JKurucz fc AvrettHigSlh . 

As regards the effective temperatu re and gravity we 
confirm the values found by ICatanzard ([2008) (see Fig. 1 
of that paper): T^ff = 24000 ± 250 K and logg = 3.91 ± 0.05, 
while ^ = 8.1 ± 0.9 km s~^ and the abundances listed in Ta- 
ble[2l 

As a ge neral result, we deri ve abundances close to the 
solar values (JAsplund et al.ll2005l ) , with the exception of car- 
bon and iron for which a slight under-abundance has been 
inferred. The slight under-abundance of magnesium, here 
derived using only one line, has to be confirmed. 



Table 2. We report the atmosphere parameters of (3 
Cep here determined: T^ff, logg, microturbulence (g) 
and abundances. N is the number of identified lines for 
any elemen t. For compari son, we report t h e va lues by: 
Morel etal.] i200dt) (M06'). iHevnderickx et al] lll994l ') (H94), 
Niecmzura fc Daszvnsk a-Daszkiewwiczl l|2005l) (N05) and the 



solar abundances by lAsplund et al.l 1120051 ) (A05). 





This study 


M06 


H94 


N05 


Teff (K) 

log 9 

5 (km s~ 


24000±250 
3.91±0.05 
1) 8.1±0.9 


26000±1000 

3.70±0.15 

6±3 


24550 
3.772 


24150±350 
3.69 


El N 


This study 


logA^Ei/A^tot 
M06 




A05 



C 8 


-3.92 ±0.05 


-4.01 ±0.10 


-3.64 ±0.05 


N 31 


-4.20 ±0.07 


-4.12±0.13 


-4.25 ±0.06 


O 41 


-3.39 ±0.07 


-3.56 ±0.14 


-3.37 ±0.05 


Ne 2 


-4.00 ±0.13 




-4.19 ±0.06 


Mg 1^ 


-4.85 ±0.10 


-4.72 ±0.21 


-4.50 ±0.09 


Al 5 


-5.80 ±0.30 


-6.02 ±0.16 


-5.66 ±0.06 


Si 10 


-4.70 ±0.19 


-4.92 ±0.23 


-4.52 ±0.04 


S 3 


-4.90 ±0.20 


-4.89 ±0.37 


-4.89 ±0.05 


Fe 14 


-4.81±0.13 


-4.79 ±0.23 


-4.58 ±0.05 


" Mgii A 4481 A 







Table 3. Theoretical equivalent widths in mA for three Si OI 
lines, computed by means of WIDTH9, for three values of the 
microturbulence. Computations are based on the previous atmo- 
sphere parameters of /3 Cep. 

Wavelength (A) ^ = 5.4 6.9 8.4 (kms"!) 



4552 
4567 
4574 



160 
141 
104 



183 
161 
118 



206 
180 

131 



5 VARIABILITY OF ATMOSPHERE 
PROPERTIES 

The spectral synthesis of the average spectrum of /3 Cep 
shows the presence of strong oxygen lines in the wings of 
the H7 line. The equivalent width variability of oxygen lines 
with the pulsation period (Fig.|3} would complicate and false 
the result of any attempt to ascertain the variability of TcS 
and log (? on the basis of this line. For this reason, we prefer 
to base this study on the H/3 line only. 

Before determining the effective temperature and grav- 
ity of any single spectrum, we check evidence for a vari- 
ability in the H/3 profile in addition to the radial velocity 
one: each spectrum has been Doppler corrected for the ve- 
locity determined by means of silicon lines and divided by 
the average H/3 profile. These ratios are shown in the upper 
panel of Fig. [2] The usual spectral distortion due to pulsa- 
tions characterises only the line core; this is here interpreted 
as evidence for a larger velocity presented by the outer lay- 
ers. The variability of the radial velocity determined through 
the core of the H/3 line presents an amplitude of 18.7±0.3 
kms~^ and an average value of — 3.3±0.2 kms~^. 



5.1 Microturbulence 

To point out any possible variability of stellar parameters 
with the pulsational period, we applied the previous iterative 
procedure to any single spectrum of /3 Cep. 

If we consider the effective temperature, gravity and mi- 
croturbulence as free parameters, the microturbulence from 
nitrogen lines along the pulsation cycle results variable be- 
tween 5.4 and 8.4 km s~^. Differently, the oxygen lines give 
a constant value equal to 8.1±0.9 km s~^. 

Source functions of the considered 14 nitrogen and 17 
oxygen lines are equally distributed in optical depth between 
logr ° = —0.85 and —2.06. These rule out any changing 
of the microturbulence with the atmospheric height. 

The number and equivalent widths of the silicon lines in 
our spectra are not enough to measure the microturbulence. 
However, the Sim triplet at 4552, 4567 and 4574 A is very 
sensitive to this velocity field and rules out a microturbu- 
lence range larger than 1.0 km s~^ . Table|3]shows that for the 
microturbulence variation deduced from nitrogen lines, the 
silicon triplet would present a 25% variability of equivalent 
widths much larger than the 7% observed (Fig.|3}- We con- 
clude that the microturbulence variability determined from 
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4855 4860 4865 4870 

Wavelength (A) 



Figure 2. The H/3 profile showed in the bottom is the average 
of single profiles acquired through a complete oscillation period. 
The upper panel shows the ratios among single profiles and the 
average one. The 1 % error bar is also reported. 
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Figure 3. Equivalent width of representative oxygen, nitrogen 
and silicon lines are reported to show their variability with the 
pulsation period. It is worth noting that O II changes out of phase 
compared to N II and Si III lines, while it is in phase with Si IV. 



equivalent widths of the representative lines of Fig. [3] for at- 
mosphere models with 20000 sC T^s ^ 30000 K, logg = 3.91, 
^ = 8.1 km s~^ and abundance listed in Table 2. It appears 
(Fig.[T} that the observed behaviour of these lines is due to 
a variation of the effective temperature betwee n 25000 and 
26000 K. That is a result in agreement with I Morel et al.l 
(120061 ) who derived an effective temperature of 26000 K from 
the ratio of Si iv and Si ill line strength. The top panel of 
Fig. [T] shows the N ll abundance computed along the pul- 
sational phase, considering the corresponding values of TcS 
and log g and keeping constant the value of ^ (8.1 km s~^). 
This reflects the variations of the ionisation state of nitro- 
gen. 



nitrogen lines is a local minimum of the previous x^ func- 
tion and we assume a constant value for this velocity equal 
to the one we obtained from the average spectrum: ^ = 8.1 
km s~^. 



5.2 Tcff and logg variability 

We find that the effective temperature and gravity of /3 Cep 
are variable with the pulsation period (Tab.[T|. Fig.[T]shows 
that Tcff changes, almost sinusoidally from ~23700 K to 
24400 K, in phase with the logg variability from 3.76 to 
3.93. 

For a constant 5, the variability of equivalent widths 
with the pulsation period is also indicative of the temper- 
ature behaviour along a cycle. Fig.O shows that N ll and 
Si III lines are in phase and that they change out of phase 
compared to O ll and Si iv. We have computed the expected 



6 CONCLUSIONS 

With the aim of determining the variations of the effective 
temperature, surface gravity and microturbulence velocity 
that occur in hot pulsating stars during the pulsational cy- 
cle, we performed time resolved R = 20000 spectroscopic ob- 
servations of the prototype /3 Cep along a complete period 
coverage. 

As a first step, we confirm t h e Teff and logg found 
in a previous paper by ICatanzarol l|2008r ) and we measure 
^ = 8.1 ±0.9 km s~^, from a sample of single lines of On. 
With these values, the abundances resulted to be solar for 
nitrogen, oxygen, neon, aluminium, silicon and sulfur, while 
slight under-abundances have been inferred for carbon, mag- 
nesium and iron. 

We find that along a pulsation cycle, after a correction 
for the pulsational Doppler shift, the H/3 line profile changes 
are matched assuming a variability of the effective temper- 
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Figure 4. Equivalent widths as a function of effective tempera- 
ture. 



ature in the range of 23700 to 24400 K and a variability 
of gravity in the range of 3.76 to 3.93. The amplitude of 
« 0.08 dex measured in the gravity variation could lead the 
reader to the erroneous belief that it correspo nds to a ra- 
dius v ariation larger than the one computed bv lAerts et al.l 
([199J), via cinematic considerations. However, the variabil- 
ity of log g here determined for /3 Cep is not a consequence 
of changes in the surface gravity because of expansions and 
compressions of the star as a whole. We ascrib ed the logo 
variability to a change of the electron pressure (lGravlll992l ') 
induced by temperature, that modulates the ionisation state 
of metals. Thus, it should not be linked directly to the radius 
of the star. 

Metal lines are characterised by equivalent widths vari- 
able with the pulsation cycle. These are explained by a vari- 
able effective temperature between 25000 and 26000 K as a 
consequence of the changed ionisation state. It could be that 
the most extended formation region of the H/3 line is at the 
origin of the lower temperature, in agreement with the 15% 
larger amplitude presented by the radial velocity measured 
through the H/3 core. 

As to the possibility that the microturbulence depends 
on the compression and expansion of the atmospheric layers, 
we find that it is constant within a 1.0 km s~^. 

We conclude that pulsations modify the dependence of 
the temperature and pressure on the optical depth in the at- 
mosphere of /3 Cep. This mimics a variability of the effective 
temperature and gravity when a small range of wavelengths 
is selected. We hope that this information can help in mod- 
eling the variability of spec tral lines that is not only due to 
the pulsation velocity field l|Dupret et al.ll2002l ). 

Along the expansion phase, the negative radial velocity 
is expected to increase up to a null value just in coincidence 
with the maximum of the stellar radius. At this phase of 
minimum compression, the photosphere is expected at the 
temperature minimum. Differently, we observe a delay in the 
temperature minimum about one tenth of the pulsation pe- 



riod. According to lDe Ridder et al.l l|2002r ). this extra-phase 
shift is due to the non-adiabatic condition that occurs in the 
pulsating layers. 
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